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We investigate excited state contaminations in a direct computation of the nucleon d-terms. This 
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T^j- ■ is an important source of systematic effects that needs to be controlled besides the light quark 

mass dependence and lattice artefacts. We use maximally twisted mass fermions with dynamical 
light (u,d), strange and charm degrees of freedom. Employing an efficient stochastic evaluation 
, of the disconnected contribution available for twisted mass fermions, we show that the effect of 

excited states is large in particular for the strange <7-terms, where it can be as big as <9(> 40%). 
This leads to the unfortunate conclusion that even with a source-sink separation of ~ 1.5 fm and 
a good statistical accuracy it is not clear, whether excited state effects are under control for this 
£3 ■ quantity. 
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1. Introduction 

The various evidences for the existence of dark matter have led to the development of exper- 
iments dedicated to detect dark matter directly. The detection relies on the measurements of the 
recoil of atoms hit by a dark matter candidate. One popular class of dark matter models involve an 
interaction between a WIMP and a Nucleon mediated by a Higgs exchange. Therefore, the scalar 
content of the nucleon is a fundamental ingredient in the WIMP-Nucleon cross section. In this 
way, the uncertainties of the scalar content translate directly into the accuracy of the constraints on 
beyond the standard model physics. Since the coupling of the Higgs to quarks is proportional to 
the quark masses, it is important to know how large scalar matrix elements of the nucleon are, in 
particular for the strange and charm quarks. 

One common way to write the parameters entering the relevant cross section are the so-called 
sigma- terms of the nucleon: 



where m denotes the light quark mass and m s the strange quark mass. To quantify the scalar 
strangeness content of the nucleon a parameter is introduced, 



The direct computation of the above matrix elements is known to be challenging on the lattice 
for several reasons. First, it involves the computation of "singlet" or "disconnected" diagrams that 
are very noisy. Second, discretisations that break chiral symmetry generally suffer from a mixing 
under renormalization between the light and strange sector, which is difficult to treat in a fully 
non-perturbative way. 

However, as has been shown in [[!]], twisted mass fermions offer two advantages here: they 
provide both an efficient variance noise reduction for disconnected diagrams ^ and avoid the 
chirally violating contributions that are responsible for the mixing under renormalization. 

A great effort has been spent developing techniques to estimate efficiently the relevant dis- 
connected contribution (see for instance ® |], [5], ^]). Also excited state contributions to nucleon 
matrix elements connected to deep inelastic scattering received a lot of attention during the past 
years [0, [||]. However, the determination of systematic effects and in particular the excited state 
contaminations of the nucleon a-terms are so far quite limited and it is the main goal of this con- 
tribution to investigate these excited state effects. 

2. Lattice Techniques 

In this study we used a single set of gauge configuration produced by the ETM collaboration. 
We used aA^ = 2+ l + l ensemble with a pion mass of 390MeV, a volume V = 32 3 x 64 and a 
lattice spacing a « 0.078 fm. We refer to [Q] for details on the gauge ensemble used in this work. In 
order to compute matrix elements involving strange quarks, we work within a mixed action setup 



o nN = m(N\uu + dd\N) and a s = m s (N\ss\N) , 



(1.1) 



2(N\ss\N) 



(1.2) 



{N\uu + dd\N)' 

which can be also related to the sigma terms of the nucleon in eq. ( |1.1| ). 
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introducing an additional doublet of degenerate twisted mass quarks of mass jj, q in the valence 
sector. 

The scalar matrix elements involved in Eq. can then be computed using the asymptotic 
behaviour of a suitable ratio of three and two-point functions denned as 

Ro q (t s ,t op ) = ^^y^ = (N\O q \N)^ + 0{e-**) + e(e-^-'^) , (2.1) 
where O q refers to the operator in which we are interested in, namely 0\ = uu + dd and O s = ss. In 



eq. (gjj) t s refers to the source-sink separation and t op to the source-operator separation. In addition, 



A represents the mass gap between the nucleon and its first excited state. From eq. ( |2.1[ ) it is clear 
that large times t op and t s are needed to suppress the so-called excited state contributions. However, 
due to the exponential decrease of the signal-over-noise ratio at large times, it is numerically very 
expensive to obtain a good signal for increasing t s or f op . 

The nucleon states themselves are created using smeared interpolating fields, which have al- 
ready been optimized themselves to suppress excited state contaminations in the two point func- 
tion. Using the same gauge field ensemble and interpolating fields, it was shown in []|] that the 
axial coupling of the nucleon, can be safely extracted with t s = Yla ~ 0.9 fm. 



Instead of using the time dependence in Rq c of eq. (|2JJ) and looking for a plateau behaviour, 



another way to extract the desired matrix element is to consider the so-called summed ratio method. 



Integrating Eq. ( |2. 1| ) over the time of insertion of the operator we are left we the following asymp- 
totic behaviour: 

Po q (t s ) = £ R 0q (t s ,t op ) =A + (N\O q \N)^t s + ff(e- At >) . (2.2) 

'op=0 

For the precise expression of the operators O q , their multiplicative renormalization and our 
computational techniques we refer the reader to [j|]. As presented in ^ part of the computation 
has been done making intensive use of a modified version of the QUDA library [M P]. 



3. Excited state contaminations: <7-terms 

We first concentrate on the computation of (Af|.w'|A , '}( bare ). Note that once multiplied by the bare 
strange quark mass, we obtain the strange a— term of the nucleon o s defined in Eq. ( |1.1| ). Since 
(A^|^!s |A^) ( bare ) only involves disconnected diagrams, with our techniques we can freely change the 
source-sink separation t s . 

In the following we will refer to a small statistics sample (ss 800 configurations) as sample A 
while a large statistics sample (« 3000 configurations) will be referred to as sample B l . Note that 
we checked that the number of noise vector used to estimate the disconnected loops is large enough 
to ensure that the statistical errors are dominated by gauge noise. 

In Fig. we show the ratio Ro, obtained at a fixed value of t s as a function of ? op . We show 
the results for t s = Yla on the sample A and for t s = 18a on the sample B. In the graph we show 
constant fits of the ratio R 0s as function of t op choosing a fitting interval ? op E [t s /2 — l,t s /2 + 1] 

'The exact number of configurations used is 853 and 2921, respectively. 
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Figure 1: The ratios Ro s (t s = 12a, t op ) and Ro s {h = Figure 2: Values for the plateaux in t op of eq (2.1) 
18a, f op ) measured with a statistics of w 800 and ps as a function of t s , for the two sets of measurements 
3000 configurations as a function of f p for fixed t s . A and B. The points have been shifted vertically for 

better readability for the sample B. 



(fop G + \) — 2, (f s /2 — |) + 2] in case t s is odd). Comparing samples A and B, we observe 
an increase of the so fitted plateaux values of ~ 40%. This indicates that source-sink separations 
of ~ 0.9 fm are not large enough to have negligible systematic errors from excited states. It is 
important to note that although we have chosen a quite large source-sink separation for sample 
B (t s = 18c?) we reached a statistics giving comparable statistical error to the one of sample A 
(t s = ISa). In order to obtain a better understanding of the systematic errors coming from the 
excited state contaminations, we analyze the strange quark content of the nucleon using both the 
ratio and and the summed ratio method. 

As mentioned above, our fit interval in ? op for R 0s is always chosen symmetrically around t s /2. 
Therefore, for each choice of t s we obtain different plateaux values providing in this way a function 
F(t s ). Ideally, when effects of excited states are sufficiently suppressed, F(t s ) reaches a plateau 
value which would correspond to the true bare matrix element {A^l^s"! A'") ( bare ) . We show in Fig. || 
F(t s ) as a function of the source-sink separation t s for the samples A and B. For better readability we 
shifted vertically the results obtained from sample B. Looking first only at sample A, it seems that 
F(t s ) indeed reaches a plateau value for t s > 11 from which we would conclude that the asymptotic 
value has been found. However, going to sample B which has a higher statistics, it becomes clear 
that F(t s ) still increases beyond t s > lla and may only reach a plateau like behaviour for t s > 16. 
This finding provides a serious warning. In fact, we cannot exclude that even with a larger statistics 
the required value of t s would be larger to observe a clear plateau. 

Let us now look at an analysis of the strange quark content using the plateau summation 
method (see Eq. (2.2)). We show in Fig. |3]the summed plateau Po s defined in Eq. (12) as a function 
of the source-sink separation for the two samples of measurement A and B. For better readability 
we shifted Po s vertically in the case of the sample of measurements B. We recall that for large t s , 
the bare matrix element can be extracted from the slope of Po s - We also show using black dotted 
lines linear fits of Pq s choosing the fitting interval t s G [t\,tz] such that we obtain a confidence level 
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Figure 3: Summed ratio Po s as a function of the 
source-sink separation for the two samples A and B. 
The results on sample B have been shifted vertically 
for better readability. 



Figure 4: Fits of the summed ratio for a fixed value 
of ?2 = 20 as a function of t\ for the two samples A 
andfl . 



larger than 90%. The corresponding slopes are 0.42(8) and 0.60(10) respectively in the case of 
samples A and B. 

In order to study the excited state contaminations we define a function S(ti,tz) as the slope 
of a linear fit of Pq in the fitting window [fi , ^2] • By construction, for ti large enough, S should 
be constant as a function of t\ and then will give the bare matrix element (./V|5 r .s|./V)( baie ). In the 
following discussion we will fix ?2 to ?2 = 20a, which appears to be large enough. We show in 
Fig. |] S(t\ ,?2) as a function of t\ for fixed ?2 = 20a and for the two samples A and B. While the 
results obtained for sample A seem to saturate for t\ > 6, the results for sample B are constant 
only for t\ > 14. Thus we reach a very similar conclusion as discussed above for the ratio Rq of 
eq. ( |2~I| ) that a too small source-sink separation and insufficient statistics can be misleading and 
that a true plateau has not been reached yet. Even with the here employed high statistics to obtain 
a signal at large values of t s , it is not clear whether the excited state contributions are negligible. 

Thus, in order to have a convincing evidence that a source-sink separation of t s = 18a ~ 1.5 fm 
is sufficient, it would be necessary to analyze for an even larger source-sink while keeping a similar 
statistical accuracy as the one obtained here. We stress that with our choice of nucleon interpolating 
fields we do not find a similarly large contribution of excited states for other observables, for 
instance the nucleon axial coupling g&- 

In the light quark sector we observe an effect of excited state contamination of the same order 
of magnitude in the corresponding disconnected part. However, here the disconnected part con- 
tributes only about 10% of the bare matrix element. In order to check the size of excited state 
contributions in the connected part, we have analyzed the matrix element for two source-sink sep- 
arations, t s = 12a and t s = 16a. We show in Fig. |5| the results for the bare matrix element for 
these two values of the source-sink separation. We observe a change of the plateau value by about 
~ 10%. Note that even if the effect of excited states is thus smaller than in the strange sector, it is 
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still large compared to our statistical errors. 



4. Excited states contamination : yu parameter 




A t, /a=1 6 
=^ Fit result for t s /a-1 2 



Figure 5: The ratio Ro,{ts = 12a, f op ) and Ro,{ts = 
16a, f p) for the light quark content of the nucleon on 
samples A and B. 



Figure 6: Plateau value for the strangeness content 
of the nucleon for two different source-sink separa- 
tion. The large source-sink separation results have 
been obtained using a larger statistics in the light sec- 
tor (connected and disconnected contributions) and in 
the strange sector. 



Finally, we have also analyzed the contribution of excited states in the determination of the 
y>N parameter. This quantity can be obtained directly by computing the ratio of two three-point 
correlators, see [jl]|. As shown in Fig. |6|, the y^-parameter obtained for a source-sink separation of 
t s = 16a on the sample B is compatible with the fit result obtained in [jl|] at a source-sink separation 
of 12a with the sample A (represented by a grey band). Thus, in the ratio corresponding to the 
jAT-parameter a cancellation of systematic effects seems to occur and the total systematic error is at 
most of the order of the statistical error, namely ~ 20%. 



5. Conclusion 

As shown in a recent paper [|lj] twisted mass fermions at maximal twist offer an efficient setup 
to compute the three point correlators relevant for the extraction of the nucleon sigma-terms. In 
this proceedings contribution, we have performed a dedicated study of excited state contaminations 
on a single gauge field ensemble corresponding to a pion mass of approximately « 390MeV. We 
have shown that the strange a s term is particularly sensitive to excited state contaminations that 
can contribute an about 40% systematic uncertainty. Since the problem is not seen in other nucleon 
observables, this large effect cannot originate from our choice of the nucleon interpolating fields. 
Thus, we suspect that our finding is true in general, also for other formulations of lattice QCD. 
We therefore conclude that studies of excited state contaminations are essential to obtain reliable 
results. As a consequence of our study, we are presently only able to provide a lower bound of the 
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systematic error in the direct determination of the strange a s term of size > 40%. In the light quark 
sector, the systematic effects is of the the order of > 10%. We observe a cancellation of the effects 
of excited state contributions in the ratio that determines the y^-parameter which lead to a safer 
determination of this quantity. 

In our analysis, we have used two approaches, the standard plateau method and the summed 
ratio method. Independently of these two analysis methods, we find it to be necessary to use 
approximately constant statistical errors while the source-sink separation is increased to detect the 
real size of the excited state contributions. Such a strategy has, of course, the drawback to increase 
exponentially the number of measurements needed, but it seems to us to be a mandatory step to 
provide reliable values of the scalar quark content of the nucleon. 
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